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Abstract 

A new measurement of the mean lifetime of the r lepton is presented. Three 
different analysis methods are applied to a sample of 90000 r pairs, collected in 
1993 and 1994 with the ALEPH detector at LEP. The average of this measure- 
ment and those previously published by ALEPH is t t = 290.1 ± 1.5 ± 1.1 fs. 
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1 Introduction 



In the Standard Model, the weak charged- current coupling strength is assumed to be 
the same for all fermion doublets. This hypothesis may be tested in the lepton sector 
by comparing the rates of certain decays. The most precise universality tests involving 
the t-v t doublet are presently obtained from measurements of the decays r~ — > e~vv, 
t~ — ► fi~~vv, and — ► e~v&, the sensitivity being limited by the experimental 
uncertainties on the r lifetime and branching fractions. 

Presented in this paper is an updated measurement of the r lifetime, based on the 
momentum-dependent impact parameter sum (MIPS) method, the impact parameter 
difference (IPD) method, and the decay length (DL) method. The MIPS and IPD analyses 
are applied to the e + e~ — > t + t~ events in which each r decays to a final state containing 
a single charged track ("1-1 topology"). The MIPS measurement has a small statistical 
uncertainty because the impact parameter smearing related to the size of the luminous 
region is nearly canceled in the sum of the impact parameters of the two daughter tracks. 
The result is, however, sensitive to the assumed impact parameter resolution. On the other 
hand, the IPD measurement is subject to a statistical error from the size of the luminous 
region but is insensitive to the assumed resolution. The DL method yields a precise 
lifetime measurement from r's decaying into three-prong final states. The dominant source 
of uncertainty is the statistical uncertainty related to the natural width of the exponential 
lifetime distribution; the size of the luminous region yields a negligible contribution to 
the lifetime uncertainty. Brief descriptions of the MIPS, IPD, and DL methods are given 
below; more details are available in [1-3]. The ALEPH measurement of the r lifetime is 
further supplemented by an analysis [|J based on the three-dimensional impact parameter 
sum (3DIP) method. The 3DIP method relies on kinematic constraints to reduce the r 
direction uncertainty in events containing two hadronic r decays. 

The data were collected at the LEP e + e~ collider, at centre of mass energies around 
the Z resonance. A r mass of m T = 1777.00to;27MeV/c 2 [0 is assumed throughout this 
paper. 

2 Apparatus and data sample 

The ALEPH detector is described in detail in [6,7]. The tracking system consists of 
a high- resolution silicon strip vertex detector (VDET), a cylindrical drift chamber (the 
inner tracking chamber or ITC), and a large time projection chamber (TPC). The VDET 
comprises two layers of double-sided silicon strip detectors at average radii of 6.3 and 
10.8 cm. The spatial resolution for r-0 coordinates is 12 /xm and varies between 12 
and 22 //m for z coordinates, depending on track polar angle. The angular coverage 
is |cos#| < 0.85 for the inner layer and |cos#| < 0.69 for the outer layer. The ITC has 
eight coaxial wire layers at radii from 16 to 26 cm. The TPC provides up to 21 three- 
dimensional coordinates per track at radii between 40 and 171 cm. The tracking detectors 
are contained within a superconducting solenoid, which produces an axial magnetic field 
of 1.5 T. Charged tracks measured in the VDET-ITC-TPC system are reconstructed with 
a momentum resolution of Ap/p = 6 x 10 _4 pt©0.005 (pt inGeV/c). An impact parameter 
resolution of 28 /xm in the r-<f> plane is achieved for muons from Z — * having at least 
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one VDET r-0 hit. 

Surrounding the TPC is an electromagnetic calorimeter (ECAL), a lead/wire- chamber 
sandwich operated in proportional mode. The calorimeter is read out via projective 
towers subtending typically 0.9° x 0.9° in solid angle which sum the deposited energy 
in three sections in depth. Beyond the ECAL lies the solenoid, followed by a hadron 
calorimeter (HCAL), which uses the iron return yoke as absorber and has an average 
depth of 1.50 m. Hadronic showers are sampled by 23 planes of streamer tubes, providing 
a digital hit pattern and inducing an analog signal on pads arranged in projective towers. 
The HCAL is used in combination with two layers of muon chambers outside the magnet 
for fi identification. 

The analysis is based on data samples collected in 1993 and 1994, corresponding to 
an estimated 31900 and 82700 produced r pairs, respectively. The combined sample of 
e + e~ -> r+r~ events consists of 3.5% collected at = 89.4 GeV, 91.4% at 91.2 GeV (Z 
"peak"), and 5.1% at 93.0 GeV. All of the off-peak data were obtained in 1993. The data 
from the two years are analyzed separately. No inconsistencies between the samples are 
observed. The results presented herein refer to the combination of the two data sets. 

Tau-pair candidate events are selected by means of an improved version of the 
algorithm described in ||. The modifications make the selection less sensitive to low- 
energy clusters in the calorimeters. The overall efficiency of this algorithm is 78%, with 
an expected background contamination of 1.3% at the Z peak. A total of 90408 candidate 
t + t~ events are selected. Further selection criteria are then imposed to isolate well- 
reconstructed one-prong and three-prong r decays for the different lifetime analyses. 

Monte Carlo t + t~ events [f| are used to study biases in the analysis methods. An 
independent sample is generated for each year of data taking to simulate the applicable 
detector conditions, beam sizes, and centre of mass energies. The input r lifetime is 
296 fs. Background events from e + e _ — > e + e~ [fT0"| , e + e~ — > fi + fi~ 0, e + e~ — > qq p"T 



and two-photon interactions |12| are also simulated. 

In the following, d denotes the impact parameter of a reconstructed daughter track, 
measured in the r-0 plane with respect to the nominal interaction point. The sign of d 
is defined to be that of the z component of the particle's angular momentum about this 
point. 

Reconstructed track impact parameters are corrected for systematic offsets in both d 
and z due to detector alignment and drift field parametrization errors. The offsets are 
measured as a function of 9 and <fi from Z — > qq events. The corrections have an rms of 
14 fim in d and 20 fim in z. Similar results are obtained from Z — > fi + fi~ events, confirming 
that the offsets do not depend appreciably on track momentum. These corrections have 
a small effect and are taken into account in the evaluation of the systematic uncertainty 
in each analysis. 

The beam axis position in the xy plane is determined from selected reconstructed 
charged tracks in Z decay events (excluding r pairs), averaged over blocks of roughly 75 
consecutive events. The typical uncertainties are 20 /jm in x (horizontally) and 10 /im in 
y (vertically). The size of the luminous region along x (typically 159 /im rms in 1993 and 
125 fim in 1994, due to different LEP optics) is measured from the fitted primary vertices 
of selected Z — > qq events, over blocks of about 270 events. The size in the y direction is 
taken to be 5 fim rms. The measured z dimensions, 7.4 mm in 1993 and 7.0 mm in 1994, 
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are used in the DL analysis. 



3 Momentum-dependent impact parameter sum 
analysis 

The MIPS method p| is applied to r + r _ events of 1-1 topology. The mean r lifetime 
is measured by means of a maximum likelihood fit to the distribution of the impact 
parameter sum 5 = d + + gL, where d + and gL denote the impact parameters of the two 
charged daughter tracks. The mean lifetime determines the overall scale of the true 5 
distribution; typical 5 values are on the order of cr T ~ 88 fim. The likelihood function 
is constructed to take into account, on an event-by-event basis, the dependence of the S 
distribution on the momenta p± of the daughters. A large sample of Monte Carlo events is 
used to parametrize this dependence for each event helicity combination, t^t^ and t^t£ . 
The likelihood function for each event includes contributions from both event helicities 
according to the measured r polarization |T3| as a function of the r~ polar angle and ^/s. 



The fitted lifetime is sensitive to the d resolution assumed in the likelihood function. 
The d resolution function for each daughter track is taken to be the sum of three Gaussian 
functions whose widths and amplitudes are parametrized || in terms of momentum, polar 
angle, and the d uncertainty calculated by the track fitting program. The parameters 
of the resolution core and near tails are obtained from studies of e + e _ — > e + e~, 
e + e~ — > yU + /i~, 77 — > e + e~, and 77 — > /i + n~ events in the data. The parameters of the 
far tails, dominated by photon conversions (which can cause pattern recognition errors) 
and nuclear interactions in hadronic r decays, are extracted from a sample of simulated 
r + r - events. The rms of the core is typically 25 to 100 /im. The second (third) Gaussian 
function contains roughly 5% (0.2%) of the tracks and is 3.3 (~20) times wider than the 
core. Finally, the small smearing in 5 related to the size of the luminous region (typically 
5 /im rms) is taken into account in the likelihood function. 

The event selection [^] yields 40271 candidate events. A few remaining events with 
poorly measured tracks are removed by means of a cut on the event confidence level (CL). 
The CL of an event is defined to be the integrated probability density for the event to 
have a reconstructed \S\ equal to or larger than the observed value; a mean r lifetime 
of 293.7 fs |§ is assumed in this calculation. The requirement CL > 0.01% removes 29 
events from the sample. The fits to the data yield an (uncorrected) mean r lifetime of 
289.1 ± 1.8 fs. Figure p] shows the 5 distribution for the data. 

The same procedure, including the parametrization of the d resolution, is followed 
with Monte Carlo e + e~ — > and 77 — > events {t = e, fi, r) generated at 

= 91.25 GeV. The different dimensions of the luminous region in 1993 and 1994 are 
taken into account in the fits to the Monte Carlo events. The possible lifetime bias in 
the method is determined from a comparison of the fitted lifetime with the generated 
value; the resulting bias is (+0.24 ± 0.63)%, where the uncertainty is from Monte Carlo 
statistics. The systematic error associated with the use of peak Monte Carlo events to 
simulate off-peak data is negligible. 

The calculated biases are used to correct the results of the fits to the data; the Monte 
Carlo fit uncertainty is therefore treated as a systematic uncertainty on the measured 
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-0.2 -0.1 0.1 0.2 

Impact parameter sum 8 (cm) 

Figure 1: Impact parameter sum distribution for 1993 and 1994 data. The curve represents 
the sum of the functions obtained from separate fits to the two data samples. 
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Table 1: Systematic biases and uncertainties in the MIPS analysis. 



Source Bias and uncertainty (%) 

t+t- Monte Carlo bias +0.24 ± 0.63 

d resolution parametrization ± 0.95 

Variation of CL cut ± 0.25 

Transverse r polarization correlation +0.22 ± 0.22 

Branching fractions —0.06 ± 0.13 

Backgrounds —0.36 ± 0.09 

Total +0.04 + 1.20 



lifetime. The systematic biases and uncertainties are summarized in Table |]. 

The systematic uncertainty associated with the d resolution parametrization includes 
contributions from (1) the limitations of the resolution model, estimated from a 
comparison of the parameters derived from simulated r + r~ vs. e + e~ and events; (2) 

the statistical uncertainty on the parameters measured from real e + e~ and n + fi~ events; 
(3) the statistical uncertainty on the parameters of the far tails, measured from simulated 
t + t~ events; and (4) the simulation accuracy of the far tail parameters, estimated from 
a comparison of real and simulated e + e~, and qq events. An additional test of the 

resolution in the r + r~ sample is performed in which the CL cut value is varied between 
and 0.5%. The resulting variations of the fitted lifetime in data and Monte Carlo agree 
within the expected statistical fluctuations; the assigned systematic uncertainty, 0.25%, 
reflects the sensitivity of the test. Any remaining detector alignment errors are implicitly 
taken into account in the d resolution parametrization. 

The correlation of the transverse polarizations of the r + and r~ is not simulated 
in the event generator || used for the final lifetime bias determination, so a special 
generator |14j] without initial or final state radiation is used to determine the effective bias 
due to this correlation, (+0.22 + 0.22)%. The uncertainty associated with the longitudinal 
r polarization is negligible. 

Each r decay mode produces a different impact parameter distribution. The differences 
between the measured r branching fractions [15,16] and those used in the Monte Carlo 
simulation are expected to yield a bias of —0.06% on the measured lifetime. The 
uncertainties on the measured branching fractions correspond to a lifetime uncertainty of 
0.13%. 

The bias due to background events is predicted from Monte Carlo simulation to 
be (—0.36 ± 0.09)%, where the uncertainty includes a systematic contribution of 25%, 
estimated from a comparison of the real and simulated distributions of the discriminating 
variables used in the t + t~ event selection. The reaction 77 — > £ + £~ is the dominant 
source of contamination in the 1-1 sample, amounting to 0.27%. The contamination from 
cosmic rays is measured to be less than 0.01%. 

The net systematic bias is (+0.04 + 1.20)%. The r lifetime result, corrected for biases, 

is 

t t = 289.0 ± 1.8 (stat) ± 3.5 (syst) fs. (1) 
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4 Impact parameter difference analysis 



The 1-1 topology events are also analyzed with the IPD method In this method the 
following quantities are determined for each event: 

Y = d + - d. , 

Ft 

where p T {\fs) is the mean r momentum, determined from Monte Carlo simulation after 
all event selection criteria are applied, p® = 45.40 GeV/c is the mean r momentum at 
y/s = 91.25 GeV, A<f) = <p + — ± it is the acoplanarity of the two daughter tracks, and 9 
is taken to be the polar angle of the event thrust axis, calculated from the reconstructed 
charged and neutral particles. No estimate of the r direction is needed to determine A0. 
For a given value of X, the expected value of Y is given by p| 



(Y) 



P°t 



X, (3) 



i.e., the slope of (Y) vs. X is equal to the mean r decay length in the laboratory frame. 
This relation holds in the approximation that the r + and r~ are back to back in the r-(f> 
projection and the r decay angles are small in the lab. 

Equation [3] is not satisfied for radiative Z — > t + t~ events; such events are rejected 
by a cut on the invariant mass of the charged daughter and the photon candidates in 
each event hemisphere. The parameters of the line (Y) = + a%X are then extracted 
by means of an unbinned least-squares fit, and the mean r lifetime is computed from the 
fitted slope a\. Event i is weighted by l/(AYi) 2 in the fit, where AYj, the expected rms 
smearing on Y it includes contributions from the estimated tracking resolution for event i, 
the size of the luminous region (a function of <f>±), and the natural spread of Yi (a function 
of Xi). The fit range |X| < 0.18 is chosen in order to reduce the effect of mismeasured 
tracks, radiative events, and background from two-photon interactions; 38120 events enter 
the fit. An iterative trimming procedure is used to remove events with fit residuals in Y 
greater in magnitude than A tr j m = 0.137 cm, most of which contain mismeasured tracks. 
This procedure removes 78 events. 

The fits to the data yield a x = +0.2218 ± 0.0023 cm and a = +0.0004 ± 0.0001cm 
(Fig. 0); the x 2 P er degree of freedom is 1.02, implying that AY^ is correctly parametrized. 
The small positive offset in is caused by bremsstrahlung and other track measurement 
errors and agrees with the value predicted from Monte Carlo events. 

Simulated events are used to study and correct for the bias in the determination of 
the mean r decay length from the Y vs. X distribution. The dependence of a\ on r T is 
influenced by the following effects: (1) The selection procedure may introduce a bias on 
the lifetime. (2) Surviving radiative events violate the assumption that the r + and r _ are 
back to back; they also cause the mean r momentum to depend slightly on |Jf |. (3) There 
is a bias associated with the approximation that the r decay angles are small in the lab 
system. (4) Tracking errors on d and introduce a positive bias on a 1; which is reduced 
by the trimming in the fit. (5) Background events affect the fitted slope. These biases 
and all systematic uncertainties are given in Table |2|. The net bias is (—0.29 ± 0.20)%, 
where the uncertainty is from Monte Carlo statistics. 



6 




-0.1 



2- (b) * - 

* * » . 

* ■ . ■ * 

»— , *-* * * *^ * 

-2 

i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

-0.3 -0.2 -0.1 0.1 0.2 0.3 

X 

Figure 2: (a) (Y) vs. X for the 1993 and 1994 data. The solid line represents the fit results. 
The dashed curve shows the shape of (Y) vs. X for simulated r + r _ and background events, 
(b) Plot of pulls (deviation from fitted line divided by uncertainty). 
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Table 2: Systematic biases and uncertainties in the IPD analysis. The first uncertainty is 
from Monte Carlo statistics and the second is systematic. 



Source Bias and uncertainty (%) 

Selection bias -0.02 ±0.11 

Radiative events -0.47 ± 0.03 ± 0.09 

Small decay angle approximation —0.12 ± 0.01 

Tracking resolution and trimming +0.54 ± 0.16 ± 0.35 

Backgrounds -0.22 ± 0.04 ± 0.06 

Detector alignment ±0.15 

Total -0.29 ± 0.20 ± 0.40 



In addition to the statistical errors from the Monte Carlo simulation, the following 
systematic errors are considered. The simulation of final state radiation in selected 
Z — > t + t~ decays, relevant to Bias 2 above, is verified in an analysis of isolated photons 
in data and Monte Carlo. The lifetime uncertainty associated with the simulation 
of the tracking resolution, including multiple scattering, bremsstrahlung, and nuclear 
interactions, is evaluated from studies of the r + r~ events as well as Bhabha, dimuon, and 
two-photon events in data and Monte Carlo. A systematic uncertainty of 25% is assigned 
on the bias from each background source. The uncertainty associated with detector 
alignment errors is taken to be half the effect of the d offset corrections. The quadratic 
sum of these systematic uncertainties is 0.40%. The beam position and size contribute a 
negligible systematic uncertainty. As a check of the procedure, the A trim value and the fit 
range in X are varied; the resulting variations in the fitted slope are consistent with those 
observed for Monte Carlo events and no additional systematic uncertainty is assumed. 
The total systematic uncertainty, including the contribution from Monte Carlo statistics, 
is 0.45%. 

The statistical uncertainty on r r is multiplied by 1.03 in order to take into account the 
small dependence of the trimming bias on r T . The uncertainty on is negligible. The 
corrected a\ value corresponds to 

r T = 290.4 ± 3.2 (stat) ± 1.3 (syst) fs. (4) 



5 Decay length analysis 

The DL method p] is used to measure the mean lifetime of r's decaying into three 
charged tracks. Three-prong hemispheres with Eg = ±1 are selected from the basic 
t + t~ sample. The event sphericity axis is calculated from the reconstructed charged and 
neutral particles for each event containing a candidate decay. The three charged tracks 
are required to point within 18° of this axis. Neutral particles outside of this cone are 
discarded and the sphericity axis is recalculated. This procedure avoids the large error 
on the r direction (and consequent lifetime bias) which can occur in radiative events. 

Decays with an identified electron are rejected to reduce contamination from photon 
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conversions. The decay vertex is fitted from the full three-dimensional track information 
provided by the detector. A candidate is retained only if the vertex fit gives a x 2 CL 
greater than 3%. This cut rejects 38% of the candidates in data and 31% in Monte Carlo. 
Most of the rejected decays contain a VDET hit that is assigned to the wrong track. 

For each candidate decay, the r flight distance is evaluated by means of a least-squares 
fit in which the r production and decay points are free to vary. The position and size of 
the luminous region and the position and uncertainty of the fitted decay vertex enter the 
fit; the r flight direction is constrained to be parallel, within an uncertainty of typically 
15mrad, to the event sphericity axis, calculated as described above. The r direction 
uncertainty is parametrized from simulated events as a function of the event sphericity; 
events with low sphericity tend to have low neutrino momenta and hence a smaller r 
direction uncertainty. The uncertainty on the fitted decay length is required to be less 
than 0.3 cm, whereas the typical decay length resolution is 0.06 cm. Decay candidates 
with a fitted decay length greater than 3 cm are discarded. This requirement removes 
two candidates whose reconstructed vertices coincide with the beam pipe (radius 5.4 cm). 
Figure |] shows the decay length distribution for the remaining 10076 candidates. 

The mean decay length (£) is extracted from the decay length distribution by means 
of a maximum likelihood fit. The probability function is taken to be the convolution of a 
decreasing exponential with a Gaussian resolution function. The slope of the exponential 
is adjusted for each event according to the LEP energy, such that the fitted (£) corresponds 
to s/s = 91.25 GeV. The decay length uncertainties, calculated event by event, are 
multiplied by a global scaling factor k which is free to vary in the fit. The fits to the data 
yield (£) = 0.2174 ± 0.0023 cm and k = 1.20 ± 0.02. 

The bias on the mean decay length, calculated from r + r~ Monte Carlo events, is 
(—0.81 ± 0.38)%, where the uncertainty is from Monte Carlo statistics. 

A Monte Carlo study shows that background events from Z —>■ qq have a mean decay 
length consistent with zero and yield a bias of (—0.26 ± 0.09)%; this uncertainty reflects 
the Monte Carlo statistics and a 25% systematic contribution. The contamination from 
one-prong r decays with converted photons is less than 0.01%. 

The effects of pattern recognition errors are studied by varying the vertex \ 2 CL cut 
between 1% and 5%; a systematic uncertainty of 0.42% is assigned, based on the observed 
variations in the mean decay length in data and Monte Carlo. If the CL cut is placed below 
1% a large increase in k is observed in data and Monte Carlo, indicating the presence of 
non-Gaussian tails in the decay length resolution. 

The use of a double- Gaussian parametrization of the decay length resolution yields a 
negligible change (<0.05%) in the fitted mean decay length. The systematic uncertainty 
associated with the parametrization of the r direction uncertainty is negligible. The 
uncertainties on the detector alignment parameters correspond to a negligible error on 
(£). The beam position and size contribute a negligible systematic uncertainty. 

The systematic biases and uncertainties are listed in Table [| The total bias is 
calculated to be (—1.07 ± 0.57)%. A correction for this bias is applied to the fitted 
mean decay length. The mean momentum of selected r's in Monte Carlo events at 
y/s = 91.25 GeV is 45.24 GeV/c; this value is used to convert the mean decay length 
to a mean proper lifetime: 

r T = 287.9 ± 3.1 (stat) ± 1.6 (syst) fs. (5) 
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Figure 3: Decay length distribution for 1993 and 1994 data. The curve represents the 
sum of the functions obtained from separate fits to the two data samples. 
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Table 3: Systematic biases and uncertainties in the DL analysis. 



Source 


Bias and uncertainty (%) 


t + t~ Monte Carlo bias 


-0.81 ±0.38 


Backgrounds 


-0.26 ±0.09 


Pattern recognition errors 


±0.42 


Total bias and uncertainty 


-1.07 ±0.57 



6 Conclusions 



The procedure of |17| is used to determine the optimum weights for averaging the measured 
lifetimes from the MIPS, IPD, and DL analyses. Correlations among the statistical and 
systematic errors are taken into account. The correlation coefficient of the MIPS and IPD 
statistical errors is calculated to be 0.44 ± 0.04 in Monte Carlo events. The combined 
result for the 1993 and 1994 data is 

r T = 289.0 ± 1.9 (stat) ± 1.3 (syst) fs, (6) 

with x 2 — 0.29 for 2 degrees of freedom (CL = 87%). 

The three new results are combined with the previously published ALEPH measure- 
ments, including those obtained with the MIPS, IPD, DL, and other methods from data 
samples collected in 1989-1992 [1-3] and with the 3DIP method from the 1992-1994 
data [[|. The statistical correlations involving 3DIP are detailed in ||. The combined 
ALEPH result is 

r r = 290.1 ± 1.5 (stat) ± 1.1 (syst) fs, (7) 

with x 2 — 9.1 for 15 degrees of freedom (CL = 87%). This result, the most precise 
measurement of the mean r lifetime, is consistent with other recent measurements W^\ - 

The ALEPH measurements of the r lifetime and branching fractions may be used to 



test lepton universality. For B(r -> evv) = (17.79 ± 0.12 ± 0.06)% pf, B(r -> \lvv) 



17.31±0.11±0.05)% [ITH], and other quantities from H, the ratios of the effective coupling 



constants 19 are 



9t 

and 

9t 

9e 



1.0004 ± 0.0032 ± 0.0038 ± 0.0005 



1.0007 ±0.0032 ±0.0035 ±0.0005, (9) 



where the first uncertainty is from the r lifetime, the second is from the r leptonic 
branching fraction (B(t — > evv) in Eq. |8| and B{t — ► fivv) in Eq. |9|), and the third 
is from the r mass. The measured ratios are consistent with the hypothesis of lepton 
universality. 
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